Spin quadrupoletronics: moving spin anisotropy around 
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We show that spin anisotropy can be transferred to an isotropic system by transport of spin 
quadrupole moment. We derive the quadrupole moment current and continuity equation and study 
a high-spin valve structure consisting of two ferromagnets coupled to a quantum dot probing an im- 
purity spin. The quadrupole back-action on their coupled spin results in spin torques and anisotropic 
spin relaxation which do not follow from standard spin current considerations. We demonstrate the 
detection of the impurity spin by charge transport and its manipulation by electric fields. 
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The field of spintronics is driven by the desire to use the 
intrinsic dipole moment of the electron, resulting from 
its spin, as an information carrier. The investigation 
and design of spintronic devices has made great progress 
in understanding the accumulation of spin dipole mo- 
ments, their manipulation by, e.g., current induced spin 
torques, and their readout by electrical transport mea- 
surements. Recently these studies have been extended to 
molecular scale quantum dot (QD) devices in partic- 
ular the predicted interplay of virtual tunneling, spin 
polarization and Coulomb interaction has been experi- 
mentally demonstrated in spin valve structures 0] . Fur- 
thermore, the importance of intrinsic spin anisotropy, in- 
duced by strong spin-orbit interaction, has been demon- 
strated in several measurements of transport through sin- 
gle magnetic molecules 0] and its various effects have 
been studied theoretically Q . Using various experi- 
mental techniques, it was shown that the intrinsic mag- 
netic anisotropy can even be controlled by atomic STM 
manipulation [3] , mechanical straining Q and reversible 
charging of the molecule controlled by a gate voltage Q . 
The spin anisotropy is of central importance to molec- 
ular scale spin manipulation, as it can provide an en- 
ergy barrier preventing unwanted spin reversal. This has 
been a key motivation in the field of single-molecule mag- 
netism and proposals for quantum computing with mag- 
netic molecules (icj | also rely on spin anisotropy. 

Spin anisotropy of a quantum state can be quantified 
by the average of the quadrupole moment tensor operator 



(1) 



where Si is the i = x, y, z component of the spin operator 
S. This operator is non zero only for spin values > 1 and 
traceless. Rs components appear in the spin Hamiltoni- 
ans describing the intrinsic spin anisotropy of magnetic 
molecules. In this Letter we show that a nonequilibrium 
spin- quadrupole moment (SQM) can be induced in a QD, 



which by itself is spin isotropic, by connecting it to fer- 
romagnets. We show that the SQM as a transport quan- 
tity obeys a continuity equation and locally affects the 
QD charge and spin and thereby the measurable charge 
current. Spin anisotropy stored elsewhere can thus be 
transferred to a place of interest by quantum electron 
transport. Similar to spin splittings induced by tunnel- 
ing |3|, the effect of this anisotropy can be large com- 
pared to the intrinsic one when coupling strongly to the 
ferromagnets. Moreover, the accumulated SQM can be 
controlled by electrical fields. Accumulation of a local 
SQM requires a high-spin quantum dot with 5 = 1 (or 
larger) in one of its charge states and a finite transport 
voltage. We illustrate the importance of SQM currents 
for the charge and spin transport across an archetypal 
model of such a QD "high-spin valve" structure: two fer- 
romagnets coupled to a single orbital level side-coupled to 
an impurity spin 1/2. In the nonequilibrium stationary 
state this QD system has a finite SQM in the high-spin 
charge state, in addition to noncollinear spins in the two 
successive charge states. Throughout the paper we set 
e = h = fee = 1 and distinguish between an operator A 
and its average A = {A) only when needed. 

High-spin valve model. In QDs close to resonance, two 
adjacent discrete charge states dominate the transport. 
The simplest model Hamiltonian of a QD with nonzero 
spin in both these charge states [11} is 



H = eN + UN t N i - Js-s'. 



(2) 



Here e denotes the energy of an orbital with occupation 
operators N = J2 a N a , N a = d\d a with a =t,4- an d 
spin s = | Ylcrcr' d]jCr acT i d a t with Pauli-matrix vector a. 
The electron on the orbital is side-coupled to a spin 1/2 
s' with isotropic ferromagnetic exchange J > 0. This 
may represent, e.g., a QD coupled to a magnetic im- 
purity 12 1 or a fullerene [HI], an asymmetrically gated 
double QD 14[, or even a hyperfine coupled single nu- 
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clear spin in a molecule [15(. Notably, both sign and 
magnitude of J can be tuned in nanojunctions [J, [Hj]. 
We consider the limit U > J > V b ,T > T r where 14 
is the bias voltage and T r the tunnel rate. Then the 
electron number N is restricted by Coulomb blockade 
to or 1 and the N = 1 singlet state can be neglected 
due to the strong exchange. Keeping only the ground 
states of the coupled orbital - impurity system with spin 
1/2 (N = 0) and 1 (N = 1), respectively and energy 
difference e = e — J/4, we obtain the simplest realiza- 
tion of an isotropic high-spin QD. The orbital is tunnel 
coupled to noncollinearly polarized ferromagnets (FM) , 
Hp = ^2 rkT t r kr c\ kT c rkT , with a constant, spin polar- 
ized density of the states k (DOS) ^ fc 8{e rkT — w) ~ v rT 
where r =t, i refers to the spin of the electrons quantized 
along the polarization axis n r of the respective electrode 
r = L,R. We let the length of the vector n r denote 
the relative polarization of the density of states of the 
FM: |n r | = (y^ - v r ±)/(v r ^ + v r ±) where v r ^ > v r ± 
by the definition of n r . The tunnel coupling to the 
FMs is accounted for by H T = J2 r kar ^o-t d\c rkT + h.c. 
where a refers to the spin of electrons quantized along 
an axis fixed to the QD (the choice of which drops out 
of the calculation). The tunneling through junction r 
is assumed to conserve spin and occurs with a spin- 
independent amplitude t r , thereby setting the tunnel 
rates T r = 2ttJ2 t u rT t^.. The spin dependence of tunnel 
amplitude in Ht, t raT = {o\er l ^ Xr ' a \T)t r , arises because 
we use a different quantization axis in each part of the 
system. Here Xr = \Xr I is the angle of rotation about the 
vector Xr, which maps the QD z-axis onto the polariza- 
tion vector n r . The electrodes are held at temperature T 
and the transport is controlled (i) by biasing the electro- 
chemical potentials fi r = ±14/2 of the electrodes with 
Vb, (h) by controlling the level position e = — V g through 
the gate voltage V g , and (iii) by adjusting the relative 
polarization angle 9 (n^ • = 1 1 | cos 9). 

Transport quantities and continuity equations. The 
theory wc now develop can address the important ques- 
tion how the impurity spin in the above model can be de- 
tected by charge transport and controlled by the applied 
voltages. To understand how charge, spin and SQM can 
accumulate on the high-spin QD, we first derive the asso- 
ciated current operators and continuity equations. The 
change of the number operator TV for electrons localized 
on the QD is induced by the injected electron particle 
currents, F N = -i[H T ,N r ] where N r = Y,kr c lkT C rkr is 



the electron number on FM r: 



N 



(3) 



This follows from the conservation of the total charge of 
the system, iV tot = N + ^ r N r , in a tunneling process, 
[Ht, N tot ] = 0, and the conservation of charge on the 
QD without tunneling, [H, N] = 0. A similar considera- 
tion for the spin operators shows that the change in the 




FIG. 1. Quantum dot high-spin valve: noncollinear ferromag- 
nets coupled to an orbital level probing a spin impurity. 



QD spin S is generated entirely by the spin-conserving 
tunneling if the QD is spin isotropic, i.e. [H, S] = 0: 



(4) 



with spin-current operators Jg = — i[ifr,S r ] and spin 
polarization S r = \ ^ kTT , c\ kT a TT ic rkT ' for FM r. The 
central aspect of the high-spin QD of interest here is that 
the average of the local quadrupole tensor (fl]) can become 
nonzero: the spin triplet state is then anisotropic in addi- 
tion to spin polarized. This anisotropy is induced by the 
SQM currents entering the QD from the spin-anisotropic 
environment: 



Q- 



(5) 



This SQM continuity equation is be obtained by using the 
conservation of the total SQM of the system in the tun- 
neling, which, however, involves various canceling contri- 
butions. Instead, the SQM current operator associated 
with electrode r, = |(IgS + SI S ) - {\t s ■ S)i + h.c, 
can be obtained more easily from definition (fT]), using the 
product rule for the time-derivative and the spin continu- 
ity equation (j4]). Since (I S S) ^ (I S )(S), the expectation 
value of the SQM-current in ([5]) has no simple relation 
to the spin current. To calculate the accumulation of 
charge, spin and SQM on the QD we need to microscopi- 
cally derive a theory which describes the QD noncquilib- 
rium state determined by the competition between the 
charge, spin and SQM currents discussed above. 

Kinetic equations. The time-evolution of the (reduced) 
density operator p of the QD, describing its nonequilib- 
rium state, is determined by the QD Hamiltonian H and 
a kernel in the kinetic (generalized master) equation. We 
diagrammatically calculate the time-evolution kernel to 
the first order in T r in the framework of the real-time 
transport theory To describe the QD high-spin valve in 
the single-electron tunneling (SET) regime, while consis- 
tently accounting for quantum coherence 1_7| , all degrees 
of freedom need to be considered, including the eight 
nondiagonal elements of the QD density matrix in the 
eigen basis of H. Without making additional approxi- 
mations we rewrite the resulting kinetic equations as an 
exactly equivalent set of coupled equations for quantum- 
statistical averages of a complete set of physical oper- 
ators in which the density operator can be expanded. 
For our model these are the charge occupancies p d , p l 



(P d - 



1) of the doublet and triplet state respectively 



3 



(TV = Tr Np = p r ) . the corresponding spin accumula- 
tions S d , S* (S = TrSp = S d + S*) and the triplet SQM 
Q* (Q = TrQp = Q*). Note that in a QD high-spin 
valve the spin accumulation vectors S d , S* in the two ac- 
cessible charge states differ in orientation and magnitude 
and need to be calculated separately. Deferring the fur- 
ther technical details to elsewhere, we discuss the physi- 
cal meaning of the resulting equations in the stationary 
limit: 



= p d = - 3 7 +p d + 2 7 V - 27+ • S d + 27" ■ S 
0=p t = 3 7 + / - 2 7 V + 2 7 + • S d - 2j- ■ S* 

- h + P d + |7"P* - 37 + S d + 7~S* 
+ 2Q* - 7 - 

| 7 V - |7~P* + 47 + S d - 27"S* 
-2Q*-7" 

[|(SV+7 + SV|(S d -7 + )l] 
-[i(S* 7 -+7-S«)-i(S*. 7 -)l] 
-2 7 "Q* + Q t X0-0 x Q* 



= S d = 



= S* = 



= Q' 



(0) 
(7) 
S d x0 

(8) 
S* x 
(9) 



(10) 



Here 7 ± = J2 r 7^ and 7^ = ^T r f^(e) is the rate for sin- 
gle charge tunneling in/out of the QD through junction 
r, denoting the Fermi distribution for electrons / holes by 
/±(e) = ( e ±0-M/r + 1 )-i. Similarly, bold-faced vectors 
7 ± = J^r'Yr an d 7? = l^ n r denote the corresponding 
rates of tunneling of spin, polarized along n r , through 
junction r. Finally, = J2 r Pr is an effective magnetic 
field, with contributions from each electrode r: 



D 



0, 



1 f du 

\T r Kc J - 

-D 



dco f+(u) 
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where the cut-off needs to be set to D ~ U since we 
exclude the TV = 2 charge state of the model. This 
field represents the spin splitting induced on the QD 
by coherent virtual electron tunneling processes into the 
spin-polarized electrode r and relies on the nonzero value 
of the Coulomb charging energy U Q ■ The magnitude 
\0r\ of these exchange fields is electrically tunable Q, 
with a peak at e = fj, r and logarithmic tails. The ki- 
netic equations incorporate the conservation of probabil- 



ity, p a 



P' 



1, and the tracelessness of the SQM ten- 



sor, trQ* = J2iQli = 0- The stationary charge current 
through junction r can be calculated in a similar way: 



In = 3 7 ,V - 2 7r V + 2 7 + • S d - 2-y~ ■ S*. 



(12) 



The continuity equation ([3]) is satisfied: TV = p* = — p d = 
Iff = ha the stationary limit. 
Due to the spin polarization of the electrodes, the 
charge occupancies ©-(O and the current (fT2j) couple to 
both charge-specific spins, but are not directly influenced 
by the SQM. Eqs. ©-([1]) show that these spins couple 



back to the charge occupancies, and suffer isotropic spin 
relaxation (oc — 7 + S d and oc — 7 _ S', respectively). In 
addition there is a transfer of spin polarization from one 
charge state to the other (oc 7~S* and oc 7 + S d , respec- 
tively). The next-to-last term represents a torque on the 
charge-specific spin due to the exchange field 0. 

A central result of this Letter is that the spin accu- 
mulated in each charge state also couples to the SQM 
accumulated in the triplet state through the last term 
±2Q* -7" in ©, ©. The opposite signs of these 
terms indicate that a nonzero SQM tends to make the 
spin polarizations in the two charge states noncollinear. 
The accumulation of SQM is described by the kinetic 
equation (fTO]) : the net injection (first two terms, us- 
ing dyadic notation (ab)^ = aibj) competes with the 
isotropic relaxation (third term). The last two terms, 



(Q* x P) tj = e jU Q\ k l3i and {0 x Q*), 



1 j 



corporate a torque exerted on the SQM by the exchange 
field 0. This finite SQM results in a back-action on the 
spin when substituted into the right hand side of 
Importantly, these back-action terms are comparable to 
the other terms in ([I])-©. By solving Eq. flTUJ) for Q* in 
terms of S d and S* they can be written as 



2Q'-7" 



S A x > 



(13) 



X=d,t 



The first two terms in Q13p make the spin relaxation and 
the spin transfer anisotropic through symmetric tensors 
R\ The remaining terms have the form of a spin torque 
involving new, charge-state specific exchange fields X 
which can be shown to be noncollinear with the stan- 
dard exchange field 0. These have the two- fold effect of 
modifying the existing spin torque term in Eqs. l[8 ]) -([9" ]) . 
—> 0±0 d ' t , and adding a torque which involves the spin 
from the other charge state. The latter thus represents 
a transfer of spin torque between the two charge states. 
In contrast to the standard exchange field 0, which is 
of a purely coherent origin (cf. (TTTj) '), both X and R A 
arise from a complex interplay of dissipative and coher- 
ent processes: their (lengthy) expressions contain both 
the transition rates 7 ± and the exchange field 0. We em- 
phasize that the above spin relaxations and spin torques 
in (|13|) cannot be understood as arising from spin (or 
charge) currents: for a correct description of the spin 
dynamics the transport of SQM must be accounted for. 
By its effect on the spin, the SQM also acts back on the 
charge occupancies and the measurable charge current, 
cf. Eq. (fT2|) . From an exhaustive study of all param- 
eter regimes of the model we find that in general the 
SQM significantly affects the charge transport whenever 
spin accumulation occurs. To emphasize its impact we 
note that if one neglects the SQM in Eqs. (f6])- (fT0"l) when 
the drain electrode is most strongly spin polarized, one 
obtains an unphysical, large particle current running op- 
posite to the voltage bias direction for nearly all relative 
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FIG. 2. (Color online), (a) Stability diagram dl/dVb vs. 
V g , Vb (red = positive, dark blue = negative). Parame- 
ters: T L = 2r R , = 0.2T and D = U = 500T in units of 
the temperature T, and 9 — 0.99-7T. The high polarization 
|nx,| = |nji| = 0.99 allows all effects to be clearly illustrated, 
lower values only result in a rescaling of the current plateaus. 
(b),(c): I vs. V g for Vt = ±36T normalized to the maxi- 
mal current 7o = ((SP^a) -1 + (2Tr : l)^ 1 )^ 1 achievable for 
Vb ^ and parallel polarizations (9 = 0). Black: full result, 
red: neglecting the SQM, blue: additionally neglecting the 
spin-transfer (ST) terms in Eq. (JHJ- (JHJ . (d),(e): Spin projec- 
tions ±S d ■ n_L and =fS* • n_R affecting the current (I12|l vs. V g 
corresponding to (b) and (c), respectively. 

polarization angles 9. We now return to the questions 
raised for the studied model. 

Impurity spin detection and control. When similar 
FMs (|n£,| = |iir|) with nearly anti-parallel polarizations 
(6 « 7r) are asymmetrically tunnel coupled (Tl > Tfi) to 
the QD, the presence of the side-coupled impurity spin 
can be detected by nonlinear transport measurements. 
Fig. (2Ja-c) show that a sharp, anomalous current peak 
occurs in the thermally broadened regime where usually 
the SET current through the QD switches on. The peak 
height depends nonmonotonically on Vb and its position 
depends nonlinearly on V g . The occurrence of such a 
current peak for opposite polarities of the bias and gate 
voltage (relative to the degeneracy point) indicates the 
presence of ferromagnetic coupling of the orbital to an 
impurity spin, resulting in nonzero total spin in both ac- 
cessible charge states. Indeed, a calculation which ig- 
nores the impurity spin completely, only shows the cur- 
rent peak for forward bias. Fig. [Hd)-(e) show that in 
both cases the current peak directly measures a signifi- 
cant precession of the spin accumulation on the QD. This 
strong precession in a narrow range of voltages is possible 
due the asymmetric tunneling coupling. At the bias and 
gate voltage where the effect is maximal the magnitude of 
the field fiii(e) originating form the weakly coupled FM 
is resonantly enhanced (|e — hr\ ~ T) and matches the 
magnitude of y8i(e) from the strongly coupled FM which 
is off resonance (|e — \ 3> T). For 8 ~ it the exchange 
field /3 is a sum of the two nearly anti-parallel, equal- 
length vectors (fTTj) and is therefore perpendicular to the 
spin, which accumulates in one of the QD charge states 
opposite to either or for \ll ^ hr, see Fig.[Ud),(e). 



Although this exchange field is small, it does, however, 
cause a large spin precession over an angle ~ tt due to the 
onset of Coulomb blockade which suppresses the spin re- 
laxation. When the weakly coupled ferromagnet is tuned 
off resonance with V g or Vb, \e — /xr| T, this precession 
is switched off. The correct description of this electrically 
controlled spin resonance relies on both the spin trans- 
fer terms (cf. blue curve) as well as the SQM induced 
back-action on the spin (cf. red curve) . 

Conclusion. The illustrated interplay of nonequilib- 
rium transport of charge, spin polarization and spin 
anisotropy in high-spin nanostructures may open up new 
possibilities for electrical detection and manipulation of 
spin 0, [Il[. We believe that "spin-quadrupoletronics" 
- the transport of spin anisotropy - presents interesting 
challenges and important prospects for further theoreti- 
cal and experimental work. For instance, the design of 
experimental setups which support pure SQM currents 
(not accompanied by spin currents) is of interest for cre- 
ating spin reversal barriers in spin-isotropic systems. For 
QDs with spins S > 1 the theory can be extended to 
transport of spin-multipole moments of rank 2S required 
to describe higher degrees of spin anisotropy. We ac- 
knowledge stimulating discussions with B. Sothmann, J. 
Konig, and S. Andcrgassen. and the financial support 
from NanoSci-ERA and DFG (FOR 912). 
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